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The radiative capture reaction rate from AA to H dibaryon 
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Radiative capture rates of thermal AA + EN states into H dibaryon are calculated in the novel 
imaginary time method. The H dibaryon is assumed to be a bound state of EN with spin J n = 0 + , 
isospin 1 = 0 and strangeness —2. We consider El transition to H from EN (L = 1) scattering states 
which mix with AA(L = 1). In order to calculate the transition rates, we formulate a coupled-channel 
imaginary time method by extending the one-channel formula originally proposed by Yabana and 
Funaki. The imaginary time method allows us to avoid the sum over all the excited thermal initial 
states, and thus to save computational time significantly. The transition rates are given as a function 
of temperature and the unknown binding energy of the H dibaryon, which we take as a parameter. 
It is found that the transition rate is not sensitive to the choices of the H binding energy or the 
strengths of the channel coupling for temperatures 3 MeV or higher. 


I. INTRODUCTION 

H dibaryon is a compact bound (or resonance) state 
of six quarks, u 2 d 2 s 2 , with spin= 0, isospin= 0 and 
strangeness= —2, proposed by Jaffe in 1977[1] . The 
lowest two baryon threshold with these quantum num¬ 
bers is the S-wave, spin-singlet AA (L = 0, J = 0) state. 
The reason why the H dibaryon is interesting is that the 
quark model dynamics, which describes the ground state 
hadrons very well, suggests a strong attractive force in 
this channel. In particular, the flavor SU(3) singlet con¬ 
figuration of the six quarks is favored most by the color- 
magnetic force (oc — (<r, • rr,)(A,; ■ Xj)), which is responsi¬ 
ble for the spin splittings of the pseudo-scalar and vector 
mesons as well as the octet l/2 + and the decuplet 3/2 + 
baryons. It is also shown that the Fermi-Dirac statistics 
for quarks allows all six quarks can occupy the lowest or¬ 
bit making a fully symmetric orbital state and therefore 
no Pauli blocking effect may induce short-distance repul¬ 
sion among the quarks. Thus a compact six-quark-like 
state is expected in the quark models [2-4]. 

The quantum chromodynamics (QCD) also allows 
color-singlet six-quark states. Recent lattice calculations 
have indeed suggested existence of either a bound or reso¬ 
nance state around the AA threshold, although the phys¬ 
ical quark-mass point calculations are yet to come[5, 6]. 
Against all these theoretical indications, long experimen¬ 
tal efforts in searching the H dibaryon were unsuccessful 
so far. The observation of the double hypernucleus has 
given an upper bound for the binding energy [7], while an 
enhancement is observed in the AA production around 
and above the threshold. [8] . 

Under these circumstances, it is interesting to explore 
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possibilities of finding the H dibaryon in various pro¬ 
duction mechanism [9]. In this paper, we investigate the 
radiative fusion production of H, A + A —> yH, in an 
environment of finite temperature baryon gas. Suppose 
that a heavy ion collision produces significant number of 
strange baryons in a temperature T. Then through the 
radiative fusion, the H dibaryons may be produced asso¬ 
ciated with emitted characteristic photons. The goal is 
to calculate the production rate from the finite-T A gas 
into H dibaryons. In the present calculation we assume 
that our final state is a SN bound state and the radiative 
transition is through mixing of SN (intermediate) states 
with AA states. 

Our setup here is as follows. The final state is a spin- 
singlet S-wave bound state of S and N. The radiative 
decay occurs most strongly from the P-wave SN (J = 1) 
states. We consider both the spin singlet ( 3 Pi) and the 
spin-triplet ( 3 Pi) SN states. They mix with the P-wave 
AA states, which due to the Pauli principle is always 
spin-triplet ( 3 Pi). Then overall, the transition goes from 
AA( 3 Pi) to SN( 1 Pi, 3 Pi) to 7 + H(5N, 1 So). 

In order to calculate the production rate, we ap¬ 
ply the imaginary time method (ITM) [10], in which 
imaginary time is identified with inverse temperature. 
This method allows us to replace the sum of the ini¬ 
tial scattering-state distribution at finite temperature by 
solving a Schrodingier-type differential equation with the 
imaginary-time variable. The advantage of the ITM is 
that it takes into account the contributions of all the ex¬ 
cited initial states automatically without approximation 
and explicit treatment of them, for instance, discretizing 
the continuum energy spectrum. 

This advantage was demonstrated in the study of the 
triple-alpha radiative capture process [11], which plays a 
key role in synthesizing 12 C in stars [12, 13]. In partic¬ 
ular, at low temperatures the reaction rate of this pro¬ 
cess is very difficult to calculate, since direct capture of 
the three-alpha particles through their scattering states 
becomes important, but formal scattering theory of the 
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three-charged particles is not available. This is in con¬ 
trast with the situation at high temperatures, in which 
the reaction proceeds through the so-called Hoyle state, 
the resonant 0 + state at 7.65 MeV in 12 C [13, 14]. How¬ 
ever, the ITM that does not require any boundary con¬ 
dition of the three-charged particles could overcome this 
difficulty and was found to be very powerful method. 

In this study, we extend the ITM mentioned above in 
coupled-channel calculations, in which the three-channel 
initial states, AA( 3 Pi) and SN( 1 Pi, 3 Pi), are considered. 
We show that the coupled-channel ITM reduces conputa- 
tional time significantly, since it does not require solving 
scattering initial states explicitly. The calculated transi¬ 
tion rate at high temperatures is shown to be insensitive 
to the energy positions of the H dibaryon, though the H 
dibaryon is assumed to be bound in this calculation, due 
to a constraint of the present method. 

This paper is organized as follows. 

In sect. 2, we summarize the formulation of the 
coupled-channel ITM. The hamiltonian for the initial 
states is taken from the Nijmegen potentials supple¬ 
mented by an anti-symmetric spin-orbit (ALS) force. 
The ALS is necessary to mix the spin 0 SN state with 
spin 1. The final state energy, i.e., the H dibaryon mass, 
and its wave function are computed by assuming a simple 
Gaussian potential for 3 So SN system. 

In sect. 3, we present details of our calculation, solv¬ 
ing the imaginary time differential equation step by step. 
Our main results are shown, i.e., the reaction rates for 
various temperatures, contributions of the initial chan¬ 
nels, as well as dependencies of the results on the choices 
of parameters of the model. 

In sect. 4, conclusions are given. 


II. IMAGINARY TIME METHOD 


In this paper, we consider the transition process in the 
heavy ion collisions. As mentioned in Introduction, H 
dibaryon is assumed to be ^So bound state of S and N 
as the final state of the radiative capture reaction. AA 
states are thus not assumed to change into the H dibaryon 
directly by emitting a photon, though they are thermally 
distributed as the initial states. The initial states that are 
most likely to decay into the final state are then SN( 3 Pi 
and J Pi) states, which mix with the AA ( 3 Pi) states by 
the strong interaction. This transition may be uniquely 
indentified as a photon with a few tens of MeV is emitted. 

The ordinary form of the radiative capture thermal 
reaction rate accompanying an emission of a photon be¬ 
tween the initial and final states, which are labeled by 


i,c and /, respectively, is given by 
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Here /3 = l/k^T is the inverse temperature, E 7 = Ei~Ef 
is the energy of the emitted photon, and tjjic and 0/ are 
the wave functions of the initial and final states, respec¬ 
tively, in which c specifies the channel number of the ini¬ 
tial states, i.e. c = 1,2 and 3 for AA( 3 Pi), SN( 1 Pi) and 
SN( 3 Pi) states, respectively. M \ M is the multipole tran¬ 
sition operator for y-ray emission with a multipolarity 
A. Since we now consider the reaction with total angular 
momentum J — 1, a y-ray with A = 1 is emitted, giving 
the transition to the 1 So H dibaryon. ay accounts for the 
degeneracy of the initial state, i.e. ay = 2 ■ 1 + 1, and /i 
is the reduced mass of the initial state, which is taken as 
the reduced mass of SN. 

In order to obtain a form of the ITM, it is convenient 
to use the following spectral representation of the Hamil¬ 
tonian, 
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By substituting the above equation into the ordinary 
form of the reaction rate Eq. (1), the following ITM for¬ 
mula can be obtained, 
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where H is the Hamiltonian for the three-component ini¬ 
tial states and P is the projection operator that elimi¬ 
nates any bound initial states. In the present calcula¬ 
tions, the Hamiltonian does not give any bound eigen¬ 
states, which leads to P = 1. 

The transition operator that we adopt in the present 
calculations is given by the standard (leading-order) El 
(A = 1) transition operator and a sub-leading operator, 
which flips the spin of the particle, i.e., 


Mx/j, = 
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where k is the momentum of the emitted photon, p rep¬ 
resents the charge density, and m denotes the magnetic 
moment density. The matrix elements of this operator 
for the individual initial channels are given by 

(ipf\M\=iu,\ipic=i) = 0 , (5) 

e f°° 

(^ / |M A= i M |^ ic=2 ) = -~r-j= / dr rR f (r)R ic = 2 (r), (6) 

4v 7r Jo 

ak f°° 

{ipf\M x= 1 ^ 0 = 3 ) = j ^ dr rR f (r)R ic — 3 (r), (7) 


where e is the electric charge, and Ri C (r ) and Rf(r) are 
radial wave functions of the initial and final states, re¬ 
spectively. Note that the leading-order transition oper¬ 
ator cannot change the spin, while the sub-leading term 
connects the spin 1 state to the spin 0 final state. We 
here assume no exchange current operator which may 
change AA to SN in the transition. The magnetic cou¬ 
pling g can be expressed in terms of the iso-scalar part of 
the g-factors of S and N, gi and g 2 , respectively, in the 
following way, 
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to vi 2 = v 2 i and v 23 = v 32 , and no terms appear in 
V 12 = v 2 \ and v 23 = v 32 other than the ALS term. 
Furthermore, one can prove that under SU(3) symme¬ 
try limit the following relation is satisfied, 

VALS = V 12 = V 21 = V 23 = v 32 = -2uan, (13) 

where van is the ALS potential between AN(‘ ? Pi) and 
AN( 1 Pi) states. 

Next we calculate the final-state wave function by 
solving the Schrodinger equation, assuming that the H 
dibaryon is a SN bound state. Here we simply adopt 
a Gaussian form for the potential of the H dibaryon as 
follows: 

VH = K>e- (r/rH)2 , (14) 

where rn is a width parameter corresponding to the size 
of the H dibaryon, for which we take rn = 0.5,1.0,1.5 
[fm]. The size of the width parameter rn = 0.5 [fm] 
corresponds to the compact 6 quarks H dibaryon and 
ru = 1.5 [fm] corresponds to a SN bound state with more 
dilute density structure. Figure 1 shows the relations 
between the various choices of r h and Vo, and the binding 
energies of the H dibaryon from the AA threshold. In 


Here mi and m 2 are masses of S and N particles, respec¬ 
tively, and pn is the nuclear magneton. 

The Hamiltonian H for the initial states can be de¬ 
scribed as having the three-channel components. The po¬ 
tential parts, which we denote as v cc ' with c, c' = 1, 2 ,3, 
can then be composed of the central, spin-orbit (SLS), 
antisymmetric spin-orbit (ALS), and tensor terms, as fol¬ 
lows: 


Vcc' — 
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where the SLS, ALS and tensor operators, L ■ ++ L ■ 
S and St, respectively, can be defined as, 

L-S^ =L-(S 1 + S 2 ) = L- (10) 

L ■ S {A) = L ■ (Sx - S 2 ) = L ■ , (11) 

St = 5i2 = 3 (oq • r) (er 2 • f) - cr 1 ■ er 2 . (12) 


Here we note that for the YN (hyperon-nucleon) and 
YY (hyperon-hyperon) interaction in Eq. (9) we adopt 
G-matrix-type YN and YY interaction constructed from 
Nijmegen potential. For the G-matrix-type interaction, 
which is described in terms of multi-range Gaussian form, 
we take the parameter set of NSC-97f [15] for the ALS 
term and those of ESC08c [16] for the other terms. The 
detailed list of the parameter sets is shown in Appendix. 

It should be noted that the ALS term does not ap¬ 
pear in un, v 22 , v 33 and U 13 = U 31 and only contributes 



Vo [MeV] 


FIG. 1. Binding energies for various choices of potential pa¬ 
rameters of the H dibaryon, m and Vo in Eq. (14). 


the subsequent calculations, we adopt rn = 1.0 [fm] and 
Vo = —219.44 [MeV], giving the binding energy of the H 
dibaryon Ef = 3.0 [MeV] from the AA threshold. The 
potential curve for this choice of the parameters rn and 
Vo is shown in Fig. 2. A radial grid size of Ar = 0.01 [fm] 
is taken in all the calculations in this work, for which we 
confirmed that the present results sufficiently converge. 
The radial wave function of the H dibaryon with this 
potential is also shown in Fig. 3. 

In order to calculate the reaction rate via the ITM 
in Eq. (3), it is useful to first define the follow¬ 
ing wave function with three components M7 (/3) = 

(4>i(/3), T 2 (/3), 4' 3 (/3)) depending on the inverse tempera¬ 
ture /3 that is identified with the imaginary-time variable, 

* m = 


(15) 
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method, in the following way, 
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fcmax ^ — Ash') 

-£ L ifer L * G9 >’ 


fc =0 


(18) 


where we take A/3 = 1 x 10 -6 [MeV -1 ] and £: max = 4. 
We confirmed that with this choice of the variables the 
results well converge. 


III. RESULTS AND DISCUSSION 


FIG. 2. The potential of the H dibaryon with Vo = —219.44 
[MeV] and th = 1.0 [fm]. The horizontal axis represents a 
relative distance between the H and N particles. 



FIG. 3. The radial part of the wave function of the H 
dibaryon. The horizontal axis represents a relative distance 
between the S and N particles. 


This wave function is found to satisfy the following 
Schrodinger-type equation along the imaginary-time axis 

P, 

- Ay (£) = h* (/3). (16) 


The formula in Eq. (3) can then be expressed by using 
the imaginary-time wave function \F(/3) as follows: 
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The reaction rate at an arbitrary value of /3 can be ob¬ 
tained as a result of the imaginary-time evolution of \F(/3) 
by solving Eq. (16) with the use of the Taylor expansion 


Figure 4 shows the initial distributions with /? = 0, 
of the AA, 5N( 1 Pi), and SN( 3 Pi) components of the 
\F(/3/2 = 0) wave function in Eq. (15), Ti(/3/2 = 0), 
Vp2(/3/2 = 0) and 4/3(/3/2 = 0), respectively. We can 
see that there is no AA component, i.e. 4'i(/3/2 = 0) = 
0, since the transition operator M\- i M does not couple 
the final state with the AA initial state, as calculated 
in Eq. (5). The amplitudes of SN( 1 Pi) and SN( 3 Pi) 
components are non-zero, and the former is larger than 
the latter, together with a broader shape of the former 
component, since the factor k in the transition operator 
for the former channel in Eq. (6) enhances the amplitude 
of the wave function at higher temperature region. 

Figures 5 and 6 show the behaviors of the imaginary- 
time evolution of the AA( 3 Pi) and SN( 3 Pi) compo¬ 
nents of the 1 ® r (/3/2) obtained by solving Eq. (16) with 
Eq. (18). In both figures, we find that the amplitudes of 
the wave functions are exuded toward the outer region 
as the imaginary-time evolution. While the amplitude of 
the AA( 3 Pi) wave function d/2(/?/2) becomes larger, as 
the imaginary-time evolution (see Fig. 5), the one of the 
SN( 3 Pi) wave function becomes smaller (see Fig. 6). 

Figure 7 is the calculated reaction rate in Eq. (17) 
shown as a function of the inverse temperature /3 [MeV], 
where the imaginary-time evolution of Eq. (16) and 
Eq. (18) is calculated in a large box size i? max = 20 
fm. We check the convergence of the reaction rate with 
respect to the box size i? max . In Fig. 8(a), we show the 
reaction rates calculated with the maximum size between 
Umax = 10 fm - 80 fm. We find that the reaction rate 
falls off too rapidly if the radial box size is not suffi¬ 
ciently large. This convergence behavior is the same as 
calculated in the previous articles with the use of the 
ITM [10, 11]. As the imaginary time proceeds, i.e., the 
temperature decreases, the amplitude of the wave func¬ 
tion gradually extends toward the outside region, since 
the lower energy states give significant contribution to 
the reaction rate at lower temperatures. We can see that 
the reaction rate is well converged up to /3 = 10 [MeV -1 ] 
if we take more than i? max = 60 fm. In Fig. 8(b), the con¬ 
vergence behavior up to /3 = 1 [MeV -1 ] is shown, and we 
can see that adopting i? max = 20 fm gives sufficiently 
converged result. Since in this study we suppose heavy 
ion collision with thermal AA gas with /3 much less than 
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(a)The AA component 4' c=1 (/? = 0). 
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(b)The SN( 1 Pi) component 4' c =2(/3 = 0). 



(c)The SN( 3 Pi) component of 4' c= 3(/3 = 0). 

FIG. 4. The radial part of the three-component wave function 
in Eq. (15) at /3 = 0. 


1 [MeV 1 ], we adopt R max = 20 fm in the subsequent 
calculations. 

In Fig. 7, we find that the reaction rate of the H 
dibaryon is 0.3779 [MeV-fm 3 ] at /3 = 0.01 [MeV -1 ] 
(ksT = 100 [MeV]). Now we can estimate the number 
of the H dibaryon in the heavy ion collision from this 



r [fm] 

(a)The wave function at f) = 0.0001 MeV -1 . 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

r[fm] 

(b)The wave function at /3 = 0.001 MeV -1 . 



FIG. 5. The radial part of the AA-component wave function 
«G=r(d/2). 


reaction rate. Suppose that one heavy-ion collision pro¬ 
duces 10 A particles within the spatial size of (5[fm]) 3 
for the time of 10/c [s], with c the light velocity. Then 
the created number of the H dibaryon is calculated to be 
N = 1.532 x 10 -3 . This indicates that for approximately 
650 collisions of heavy ions, one H-dybarion is produced. 
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(a)The wave function at ft = 0.0001 MeV 1 . 



(b)The wave function at /3 = 0.001 MeV 1 . 



(c)The wave function at ft = 0.01 MeV T 

FIG. 6. The radial part of the HN( 1 Pi)-component wave func¬ 
tion T c= 2(/3/2). 


The ITM formula of Eq. (3) or Eq. (17) could be ob¬ 
tained by transforming the ordinary form of the reac¬ 
tion rate Eq. (1), where scattering wave functions should 
be handled explicitly as the initial states. The present 
formula is the extended version of the original ITM for¬ 
mula in Ref. [10] to the case of coupled-channel calcula- 



p [MeV -1 ] 


FIG. 7. The reaction rate of the H dibaryon as a function of 
the inverse temperature f3. 



(a)Calculated for ft < 10 MeV 1 . 



p [MeV -1 ] 

(b)Enlarged view of (a) with ft < 1.0 MeV -1 . 

FIG. 8. The convergence behavior of the calculated reaction 
rates with respect to the choice of the different radial box 
sizes Umax ■ 
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tions. Therefore it is important to check the correctness 
of the analytical formula of the coupled-channel ITM, 
by comparing the results calculated numerically accord¬ 
ing to both the fomulae Eq. (17) and Eq. (1). In order 
to calculate the reaction rate according to the ordinary 
method, we first derive the initial states ipi C (c = 1, 2,3) 
with the use of the conjugate gradient (CG) method in 
a spatial box size, which is here taken as 20 [fm], with 
a radial grid size Ar = 0.01 fm. The total number of 
the initial states is then 2000 x 3. Figure 9 shows the 
convergence behavior of the reaction rate at /3 = 0.0001 
[MeV -1 ] for the adopted number of the initial states in 
the sum of the energy levels i in Eq. (1). We find that 
the reaction rate is converged when more than 400 initial 
states, i.e. c ~ V’^ooci are summed up in Eq. (1). 
The eigenenergy of the 400th energy eigenstate is calcu¬ 
lated to be -Ei= 4 oo = 15748 MeV. 



the number 

FIG. 9. The reaction rate at /3 = 0.0001 MeV -1 calculated 
with the ordinary method of Eq. (1). The convergence behav¬ 
ior with respect to the number of initial states taken in the 
sum of Eq. (1). 



FIG. 10. The comparison of the reaction rates calculated with 
the ITM and the ordinary method. 


In Fig. 10, we show the reaction rate calculated as 


a function of /3, in which the energy eigenstates up to 
*max = 400 are summed up, in comparison with that 
obtained with the ITM. We can see that both results co¬ 
incide with each other, indicating the correctness of the 
ITM formula derived within the coupled-channel frame¬ 
work. Here we mention an advantage of the ITM that 
it saves the computational time much more than the or¬ 
dinary method. The former takes 57 seconds while the 
latter 608 seconds, at /3 = 1 MeV -1 , i.e. 10 times faster 
than the ordinary method. 

We next discuss the dependence of the reaction rate 
on the binding energy of the H dibaryon, which can be 
artificially tuned by changing the depth parameter of the 
potential Vq in Eq. (14) (see also Fig. 1). Figure 11 shows 
the reaction rate with different binding energies of the H 
dibaryon. We see that as the binding energy increases the 
reaction rate becomes large at lower temperature regions. 
However, at higher temperatures, /3 < 0.3 MeV -1 , which 
we are interested in in this study, the reaction rate is not 
sensitive to the variation of the binding energy of the H 
dibaryon. This implies that a similar result is expected 
for a resonance H dibaryon above AA threshold, although 
in this calculation a bound state is assumed. 

i 
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£ 1X10" 7 
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0.0 0.2 0.4 0.6 0.8 1.0 

p [MeV -1 ] 

FIG. 11. Dependence of the reaction rate on the binding 
energy of the H dibaryon, which is varied by changing the 
depth parameter Vo in Eq. (14), with a fixed width-parameter 
value rn = 1.0 fm. 


We also investigate the dependence of the reaction rate 
on the choice of the ALS potential, by artificially vary¬ 
ing the strength of the ALS potential uals in Eq- (13). 
We show in Fig. 12 the reaction rates calculated with 
the various strength values with which additional fac¬ 
tors are multiplied. We find that at lower temperature 
region with /3 > 0.4 MeV -1 , the reaction rate is sensi¬ 
tive to the choice of the ALS potential, and the stronger 
ALS potential leads to the larger reaction rate, whereas 
at higher temperature region with /3 < 0.4 MeV -1 , the 
reaction rate is not sensitive to the choice of the ALS 
potential. This result is understood as follows: At low 
temperature region, the AA states are dominant and the 
mixing of the SN( 1 Pi) states with the AA states is small. 











The direct transition of the AA initial states into the H 
dibaryon is prohibited and only the ALS potential cou¬ 
ples the A A channel with the SN( 1 Pi) channel that con¬ 
tributes to this transition, resulting in the strong sensi¬ 
tivity to the ALS potential. On the contrary, at higher 
temperatures, there are sufficient number of the SN( 1 Pi) 
thermally mixed with AA states. The mixing of the SN 
states gives non-negligible contribution to the reaction 
rate, even without the ALS potential, leading to the in¬ 
sensitivity to the choice of the ALS potential. 



FIG. 12. Dependence of the reaction rate on the magnitudes 
of the ALS potential wals- 


ALS coupling potential to the SN( 1 Pi), and then to the 
H dibaryon. We also note that the wave-number factor 
k of the emitted 7 -ray in the transition operator for the 
SN( 1 Pi) channel in Eq. ( 6 ) makes c-n^Pi) larger than 
Chn( 3 Pi)- On the contrary, in lower temperature region, 
there are few SN states mixed with the AA thermal distri¬ 
bution. The dominant contribution to the reaction rate 
is then from the AA channel with a significant role of the 
ALS potential. 



FIG. 13. Contributions of the reaction rate from the AA, 
HN( 1 Pi), and EN( 3 Pi) channels, c AA ( 3 Pi)i CsNpPi) and 
c HN (3 P i), respectively, defined in Eq. (19). 


Finally, we calculate the contributions to the reac¬ 
tion rate from the each component of the wave function 
\I/(/3) = (Ti^), ^(/S), 'L 3 (/3)). In Eq. (17), the reaction 
rate can be decomposed into the following three compo¬ 
nents, 


c AA( 3 Pi) oc ^2 

C 

C HN(1PD OC ^2 

C 

ChN( 3 P!) OC ^2 

C 





where the following relation is satisfied, 

c AA( 3 Pi) + c HN(!Pi) + c EN( 3 Pi) = r (/3) • (20) 


We can understand this situation more clearly by the 
further analysis of the three-component wave function. 
In Fig. 14, we shows the probabilities of the each compo¬ 
nent defined below, 


PAA( 3 Pi) 


P-NpP!) 


P~N( 3 P!) 
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where the following relation is satisfied, 

PAA( 3 Pi) +P3N( 1 Pi) +P=N( 3 Pi) = 1- 


( 21 ) 

( 22 ) 

(23) 

(24) 


We show in Fig. 13 these contributions from the three 
channels. We can see that at higher temperatures the 
contribution of the SN( 1 Pi) channel, is the 

largest. This is because, according to the previous con¬ 
sideration, at higher temperatures a plenty of SN( 1 Pi) 
states are mixed with the thermal AA distribution, and 
the direct transition process from the SN to the H 
dibaryon prevails over the process from AA through the 


We can see that it is found that, as is consistent with 
the previous figure, at high temperatures the SN( x Pi) 
component is dominant, and at the lower temperatures 
the AA component becomes the largest. As the decrease 
of temperature, the probability of the AA component 
rapidly gets higher, and over /3 ss 0.4 MeV -1 it becomes 
the largest, whereas below this temperature the SN( 1 Pi) 
component rapidly disappears. 
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FIG. 14. Contributions of the probability from the AA, 
HN( 1 Pi), and HN( 3 Pi) channels, Paa(3Pi)i PsNpPi) and 
Psn(3p 1 )j respectively, defined in Eqs. (21)-(23). 


IV. SUMMARY AND CONCLUSION 

In summary, we have applied the imaginary time 
method to calculate the radiative fusion rates of AA into 
the H dibaryon at finite temperature. Mixings of the 
SN( 3 4 5 6 7 Pi) 2 and SN( 1 Pi) channels with the AA( 3 Pi) chan¬ 
nel are considered in the initial thermal state. The El 
transition to the SN( 1 So) bound state that represents the 
H dibaryon is calculated. The imaginary time method is 
applied so that the sum over all the excited states in the 
thermal initial state can be taken into account without 
computing the scattering wave functions. 

A representative calculation is done for a bound H 
dibaryon with the binding energy 3 MeV below the AA 
threshold for the initial temperature of 100 MeV, giving 
the transition rate, 0.38 [MeV-fm 3 ]. This corresponds to 
a rate that may produce one H dibaryon in a few hun¬ 
dred HI collisions through the El transition, assuming 
a collision produces 10 A’s in a few frn volume region of 
temperature about 100 MeV. It is found that the effect of 
the channel mixing is significant especially at high tem¬ 
perature, where the difference of the threshold energies 
between AA and SN is not important. On the other hand, 
at a lower temperature, the transition rates are sensitive 


to the mixing potential. It should be noted that the El 
transition is dominated from the SN( 1 P 1 ) channel, which 
mixes with the AA( 3 Pi) state only through the antisym¬ 
metric spin-orbit (ALS) potential. Thus the transition 
rates at low temperatures are sensitive to the strength 
of the ALS potential. We have found that the transition 
rate is larger for the deeper bound state at low temper¬ 
atures. In contrast, the reaction rates are not sensitive 
to the choices of parameters when the temperature is 
above around 50 MeV. Although our present calculation 
is applied only to the bound H state, we expect a similar 
production rate, even if the H dibaryon is a resonance 
above the AA threshold. 

Our calculation is compared to the standard method, 
which explicitly sums the initial scattering states, and it 
is found that the results are consistent with each other. 
On the other hand, the computational time in the imagi¬ 
nary time method is much shorter (about one order) than 
the ordinary method. Thus, we have confirmed that the 
imaginary time method for the multi-channel hadronic 
transitions can be well performed. 
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Appendix: Parameter sets of YY and YN interaction 

We show in Table I the parameter sets of the G-matrix- 
type YY and YN interaction, which are expressed in the 
Gaussian form in Eq. (9). For the central, SLS and ten¬ 
sor forces, we adopt the Nijmegen ESC08 potential, while 
for the ALS force the Nijmegen NSC97-f is chosen. The 
width and strength parameters of the Gaussian form are 
shown for each v cc > channel with c, d = 1,2,3. As de¬ 
scribed in the text, the ALS force appears only in v \2 
and U 23 , for which no terms contribute other than the 
ALS term. 
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TABLE I. The width and strength parameters of the YY and YN potential written in terms of the Gaussian form in Eq. (9) 
are shown. The parameter sets of ESC08c [16] and NSC97-f [15] are adopted for the central, SLS and tensor forces, and the 
ALS force, respectively. 


Vij 

k 

r k C) M 

v[ c) [MeV] 

r ( k s) [fm] 

v[ s) [MeV] 

r£ T) [fm] 

4 T) [MeV] 

r i A> [ fm l 

v[ A) [MeV] 

AA( 3 Pi)-AA( 3 Pi) 

1 

0.2274 

— 1.707 x 10 3 

0.05465 

-2.297 x 10 3 

0.07595 

-1.018 x 10 2 




2 

0.3213 

1.100 x 10 4 

0.09786 

-1.005 x 10 3 

0.1512 

2.113 x 10 2 




3 

0.4538 

-4.355 x 10 3 

0.1752 

-4.807 x 10 2 

0.3011 

-2.911 x 10 2 




4 

0.6411 

1.028 x 10 3 

0.3137 

2.178 x 10 3 

0.5997 

1.520 x 10 2 



Vll 

5 

0.9057 

-8.899 x 10 1 

0.5618 

-1.913 x 10 2 

1.194 

-4.329 x 10° 




6 



1.005 

1.239 x 10 1 






7 



1.801 

-7.992 x 10 _1 






8 



3.225 

7.613 x 10 ~ 2 





EN( 1 Pi)-EN( 1 Pi) 

1 

0.1673 

1.180 x 10 2 








2 

0.2778 

-3.910 x 10 3 







V22 

3 

0.4613 

1.021 x 10 3 








4 

0.7659 

5.406 x 10 2 








5 

1.271 

-3.311 x 10 1 







HN( 3 Pi)-EN( 3 Pi) 

1 

0.1817 

-1.449 x 10 3 

0.1141 

-2.052 x 10 3 

0.1456 

2.648 x 10 1 




2 

0.3165 

3.956 x 10 3 

0.2696 

-3.974 x 10 3 

0.2789 

-5.224 x 10 2 



V33 

3 

0.5512 

-3.088 x 10 1 

0.6366 

1.621 x 10 2 

0.5341 

4.469 x 10 2 




4 

0.9599 

1.982 x 10 1 

1.503 

8.211 x 10 ° 

1.022 

4.907 x 10 1 




5 

1.671 

-1.217 x 10 1 

3.550 

-4.802 x 10 ' 1 

1.958 

3.771 x 10° 



AA( a Pi)-HN( 3 Pi) 

1 

0.06457 

-9.471 x 10 2 

0.05748 

3.213 x 10 2 

0.2276 

-4.925 x 10 2 




2 

0.1308 

1.658 x 10 3 

0.1272 

6.710 x 10 1 

0.3308 

9.332 x 10 2 



U 13 (— v 3 i) 

3 

0.2653 

-4.896 x 10 3 

0.2817 

-7.415 x 10 3 

0.4806 

-5.395 x 10 2 




4 

0.5378 

8.913 x 10 2 

0.6236 

3.363 x 10 2 

0.6983 

9.507 x 10 1 




5 

1.090 

-8.803 x 10 1 

1.380 

3.634 x 10° 

1.014 

-9.035 x lO ” 1 



AA( 3 Pi)-EN( i Pi) 

1 







0.75 

-64.50 

U 12 (— V21 = V23 = V32) 

2 







0.40 

- 686.6 
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